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ABSTRACT 

An analysis of large-area CO J=3-2 maps from the James Clerk Maxwell Telescope for 
12 nearby spiral galaxies reveals low velocity dispersions in the molecular component 
of the interstellar medium. The three lowest luminosity galaxies show a relatively flat 
velocity dispersion as a function of radius while the remaining nine galaxies show 
a central peak with a radial fall-off within 0.2 — 0.4r25. Correcting for the average 
contribution due to the internal velocity dispersions of a population of giant molecular 
clouds, the average cloud-cloud velocity dispersion across the galactic disks is 6.1 ± 
1.0 km s _1 (standard deviation 2.9 km s" 1 ), in reasonable agreement with previous 
measurements for the Galaxy and M33. The cloud-cloud velocity dispersion derived 
from the CO data is on average two times smaller than the HI velocity dispersion 
measured in the same galaxies. The low cloud-cloud velocity dispersion implies that 
the molecular gas is the critical component determining the stability of the galactic 
disk against gravitational collapse, especially in those regions of the disk which are H2 
dominated. The cloud-cloud velocity dispersion shows a significant positive correlation 
with both the far-infrared luminosity, which traces the star formation activity, and the 
K-band absolute magnitude, which traces the total stellar mass. For three galaxies in 
the Virgo cluster, smoothing the data to a resolution of 4.5 kpc (to match the typical 
resolution of high redshift CO observations) increases the measured velocity dispersion 
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1 INTRODUCTION 

The vertical structure of the interstellar medium is de- 
termined by a delicate balance between gravity and pres- 
sure. The concentration of mass in the stellar disk drags 
the gas into a thin disk whose vertical scale height is con- 
trolled by the velocity dispersion. The velocity dispersion 
of the gas is also an important parameter for star for- 
mation laws based o n the Toomre Q criterion l|Toomrel 
1 19641 ; iKennicuttl 1 19891) . The atomic phase of the interstel- 
lar medium (IS M) has been well studie d and interpreted, 
most recently by iTamburro et al.l (|2009l . see below) . How- 
ever, the velocity dispersion of the star forming molecu- 
lar gas is much less well understood, although the avail- 
able data are consistent with a sig nificantly thinner and 
dynamically colder molecular disk ([g tark & Brand 1 19891 ; 
IWilson fc Scovillelll990l ; ICombes fc Becauaertlll99i r 

Determining the velocity dispersion of the molecular 
gas in the Galaxy is complicated by our location within the 
plane of the disk and the difficulty in removing the effects 
of streaming motions associated with spiral arms. Nearby, 
relatively face-on galaxies are better targets, provided suffi- 
cient sensitivity and spectral and spatial resolution can be 
achieved. However, an additional complicating factor in any 
analysis is the fact that the observed velocity dispersions are 
not significantly larger than the internal velocity dispersion 
of an individua l giant molecular cloud (GMC). For example, 
ISolomon et alj l|l987l ) measure internal velocity dispersions 
ranging from 1 to 8 km s _1 and a relation between internal 
velocity dispersion and cloud mass that goes as M = 2000<r*. 
Certainly the net effect of the internal velocity dispersion of 
a population of giant molecular clouds in the beam needs 
to be considered in an analysis of the cloud-cloud velocity 
dispersion of nearby, near face-on galaxies. 

For our own Galaxy, IClemensI (|l985h obtained a one- 
dimensional cloud-cloud velocity dispersion of 3.0 km s~\ 
sig nifican t ly sm a ller than the value of 7 -9 km s _1 obtained 
by IStarkl (l984h . Istark fc Brand! Jl989f) argue that the ve- 
locity dispersion measured by " IClemensI (|l985h is actually 
the internal velocity dispersion of individual clouds and ob- 
tain a cloud-cloud velocity dispersion of 7.8 ± 0.6 km s" 1 
for clouds within 3 kpc of the Sun. Although this mea- 
surement includes small-scale streaming motions, they argue 
that the true value is only 20% smaller whe n the streaming 
motions are removed Istark fc Brandl ll989). More recently 
IStark fc Lee] (|2005l . I2006T ) have re-derived the scale height of 
the molecular gas to be 35 pc for clouds less massive than 

2 x 10 5 M© and only 20 pc for clouds more massive than this 
limit. This scale height of 20 pc implies a one-dimensional 
cloud-cloud velocity dispersio n of just 4 km s for the 
more massive molecular clouds. ICombes fc Becquaertl (|l997l ) 
point out that, given the observed scale heights of HI and 
CO in the Galaxy and a velocity dispersion of 9 km s 
in the atomic gas, we would expect a cloud-cloud velocity 
dispersion of only 2.4 km s _1 in the molecular gas, a value 
which is significantly smaller than any of the measurements. 

There are relatively few measurements of the ve- 
locit y dispersion of the mo lecular gas in other galax- 



Wilson fc Scovi 11 (1990) used a combination of in- 
terferometric observations of individual giant molecular 
clouds with single dish observations of M33 to mea- 
sure a cloud-cloud velocity dispersion of 5 ± 1 km s . 



ICombes fc Becquaertl (|l997f l observed NGC 628 and NGC 
3938 using the CO J=l-0 and J=2-l lines and observed ve- 
locity dispersions of 6 km s" 1 and 8.5 km s~ , respectively. 
They used gaussian fits to th e observed line widths an d 
corrected for satur ation effects (|Garci'a-Burillo et al.lll993l ). 
IWalsh et all (|2002l ) made CO J=l-0 and J=3-2 observations 
of NGC 6946 and measured velocity dispersions from second 
moment maps of 8.9±2.2 km s _1 and 6.0± 1.6 km s -1 in the 
two l i nes. Neither of thes e two studies (|Combes fc Becquaerd 
1 19971 ; IWalsh et al]|2002l ) attempted to correct for the inter- 
nal velocity dispersion of individual GMCs, and so these 
measurements are upper limits to the cloud-cloud velocity 
dispersion. 

There have been a number of theoretical attempts to 
model the clou d -cloud velocity dispersion in the Galaxy. 
IJog fc Ostrikeil j 19881 ) propose that cloud-cloud gravita- 
tional scattering in a differentially rotating galactic disk acts 
to increase the random kinetic energy of the cloud pop- 
ulation. In this model, inelastic collisions between clouds 
act as an energy sink, resulting in an equilibrium value 
for the one-dimensional v elocity dispersion of 5-7 km s . 
iGammie. Ostriker fc Joel (|l99ll ) extended this work using 
both analytical and numerical analyses and obtain a value 
of 5 km s _1 for the t wo-dimensional vel o city d ispersion in 
the plane of the disk. iThomasson et all (|l99ll ) performed 
N body simulations of clouds and stars which include in- 
elastic cloud collisions but not close 2-body gravitational 
encounters. They find typical velocity dispersions of 3 km 
s _1 which incr ease in galaxies with s tronger spiral structure. 
More recently. [Tasker fc Tanl |2009l ) have developed a three 
dimensional model including ISM cooling to 300 K with a 
resolution of 8 pc and find typical velocity dispersions of 10 
km s -1 for a model of the Milky Way. 

Although there have been relatively few measurements 
of the velocity dispersion in the molecular gas in galaxies, 
the velo city dispersion of the atomic gas has been well stud- 
ied (see ITamburro et al.| [2009. and references there in). Most 
recently, the THINGS survey (| Walter et al.ll2008l ) has pro- 
duced high-resolution HI maps of 34 spiral and irregular 
galaxies with distances less than 11 Mpc. The typical ve- 
locity dispersion in the atomic gas at radii between r2s/2 
and r25 i s 11 ± 3 km s -1 for galaxies with inclinat ions less 
than 60° (|Lerov et al.ll2008l ) . ITamburro et""aH (|2009l ) find an 
HI line width that falls off systematically with radius which 
they link to the energy provided by supernovae linked to 
recent star formation. They also find a characteristic HI ve- 
locity dispersion of 10 ±2 km s _1 at T2s, which often marks 
the extent of significant star formation in the disk. 

In this paper, we present measurements of the velocity 
dispersion for the molecular component of the interstellar 
medium using data from the JCMT Nearby Galaxies Legacy 
Survey (NGLS) as well as a follow-up program on Hi-rich 
spiral galaxies in the Virgo Cluster. We use new, wide-area 
observations of the CO J=3-2 emission for 12 large spiral 
galaxies with inclinations < 60° and distances < 17 Mpc. We 
discuss the observations, data reduction, and analysis meth- 
o ds used in &\ A more detailed comparison with the results 
of ICombes fc Becquaertl (|l997l ) is given in Appendix A. In 
Jj3] we discuss the observed values of the velocity dispersion 
of the CO J=3-2 transition and the correction factor nec- 
essary due to the non-trivial internal velocity dispersion of 
individual giant molecular clouds. The potential need to cor- 
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rect for a non-isotropic velocity dispersion in galaxies with 
inclinations larger than about 30° is discussed in Appendix 
B. In iJH we compare our results with recent observations 
of the velocity disper sion of the atomic g as component from 
the THINGS survey | Walter et al.ll2008l ), investigate corre- 
lations of the velocity dispersion of the molecular gas with 
global galaxy properties such as mass and star formation 
rate, discuss the implications for understanding disk stabil- 
ity in galaxies, and compare our data with recent observa- 
tions of spiral galaxies at z — 1. We give our conclusions in 



2 OBSERVATIONS AND DATA PROCESSING 
2.1 JCMT CO J=3-2 data 

The CO J=3-2 observations were obtained as part of the 
JCMT Nearby Galaxies Legacy Survey (NGLSfl which is 
observin g an HI flux limited sample of 155 galaxies within 
25 Mpc ([Wilson et aj]|2009ft . The angular resolution of the 
JCMT at this frequency is 14.5", which corresponds to a 
linear resolution ranging from 0.2 to 1.2 kpc for the galax- 
ies in our sample. From the large (-D25 > 4') spiral galaxies 
observed by the NGLS, we selected 9 galaxies with inclina- 
tions smaller than 60° to study the velocity dispersion. Eight 
additional large galaxies from the NGLS with low inclina- 
tions (IC 2574, UGC 04305, NGC 3031, NGC 3351, NGC 
4450, NGC 4579, and NGC 4725) did not have sufficiently 
strong or extended detections in CO J=3-2 to be useful for 
this analysis. To this sample we added three spiral galax- 
ies in the Virgo Cluster (NGC 4303, NGC 4501, and NGC 
4535) which are not part of the NGLS survey but which have 
been observed in a similar manner. (A fourth Virgo Cluster 
galaxy, NGC 4548, did not have a strong enough detection 
to be useful here.) These galaxies were observed as part of a 
follow-up program to the NGLS (JCMT proposal M09AC05, 
PI C. Wilson) to obtain CO J=3-2 observations to complete 
the sample of Virgo spiral galaxies with HI fluxes > 6.3 Jy 
km s _1 . The data for these three galaxies were obtained 
between 2009 February 14 and 2009 May 26. 

Relevant properties of the 12 galaxies are given in Ta- 
ble [1] All 12 galaxies were observed in raster mapping mode 
to cover a rectangular area corresponding to D25/2 with a 
1 sigma sensitivity of better than 19 mK (T\) at a spec- 
tral resolution of 20 km s -1 . We us ed the 16 pixel array 
receiver HARP-B l|Buckle et all 120091 ) with the ACSIS cor- 
relator configured to have a bandwidth of 1 GHz and a reso- 
lution of 0.488 MHz (0.43 km s~ x at the frequency of the CO 
J=3-2 transition). The combination of uniquely high veloc- 
ity resolution and large mapping area is critical for accurate 
measurements of the velocity dispersion in the molecular in- 
terstellar medium. 

Details of the redu ction of the CO J=3-2 da ta are given 
in lWilson et al.l ([20091 ) and IWarren et alj (|2010l ) and so we 
discuss here only the most important processing steps and 
those steps that differed from the previous analysis. The in- 
dividual raw data files were flagged to remove data from 
any of the 16 individual receptors with bad baselines and 
then the scans were combined into a data cube using a 



http://www.jach.hawaii.edu/JCMT/surveys/ 



sinc(7ra)sinc(&7ra) kernel as the weighting function to deter- 
mine the contribution of individual receptors to each pixel in 
the final map. The pixel size in the maps is 7.276". A mask 
was created to identify line-free regions of the data cube 
and a first-order baseline was fit to those line-free regions 
and subtracted from the cube. 



W e then used the clumpfind algorithm (Williams et al 
I l994h implemented as part of the CUPILrl (|Berrv et al 



120071 ) task findclumps to identify regions with emission with 
signal-to-noise greater than 2.5 in a data cube that had been 
boxcar smoothed by 3x3 spatial pixels and 25 velocity chan- 
nels. Moment maps were created from the original data cube 
after applying the mask created by findclumps. The moment 
maps which are the focus of this paper are the moment 2 
maps, which measure the velocity dispersion, a v , for each 
pixel in the image using 



a v = y/YTi(vi -v) 2 /T,T t 

where Vi and Ti are the velocity and temperature of a given 
velocity channel and v is the intensity weighted mean ve- 
locity of that pixel. This method of c alculating the velocity 
disper sion differs from that used by ICombes fc BecquaertJ 
l|l997h . who fit gaussian profiles to the CO lines. How- 
ever, in the limit of gaussian lines with a high signal-to- 
noise ratio, the values from the moment 2 maps should 
agree with the results from fitting a gaussian directly to 
the l ine profiles. A more detailed comparison of our data 
with ICombes fc BecquaertJ l|l997h is given in Appendix A. 
Moment 2 maps are often used to determine the velocity 
dispersion in the atomic gas, although this method is most 
reliable for simple line shapes, as warped di sks and HI at 
large scale heights can distort the line profiles (|de Blok et al.l 
|200S| ). 

Because we were using a relatively low signal-to-noise 
threshold in creating our masks, the final moment 2 maps 
sometimes appeared to contain spuriously high values typi- 
cally in the outer portions of the map. These regions appear 
white in Figures [T] and [2] These values would bias the dis- 
persions upward if they were included in our averages. For 
5 galaxies in our sample, (NGC 4321, NGC 628, NGC 2403, 
NGC 3184, and NGC 4826), we applied an upper threshold 
to remove pixels with values that were higher than the high- 
est value seen in the central part of the galaxy. The thresh- 
olding values used were 15 km s _1 for NGC 628, 25 km s~ x 
for NGC 2403 and NGC 3184, 52 km s" 1 for NGC 4321, and 
85 km s _1 for NGC 4826. We also applied a threshold cut 
of 65 km s _1 to NGC 4501 but very high values persisted 
in the north-western portion of the map (see below); the ve- 
locity dispersions in Table [T] are measured in the south-east 
portion of the map only. For NGC 5055, we used a targeted 
threshold to remove the block of high values seen in the 
north-east portion of the map. 



2.2 HI data from THINGS 

For the 6 galaxies (NGC 628, NGC 2403, NGC 3184, NGC 
4736, NGC 4826, and NGC 5055) that are in common be- 



b CUPID is part of the Starlink llCurrie et all l200g|) 
software package, which is available for download from 
http://starlink.jach.hawaii.edu 
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Figure 1. CO J=3-2 velocity dispersion for six galaxies which are not members of the Virgo cluster, (a) NGC 628. Contour levels are 
2.5, 5, 10 km s — 1 and the colour scale peak is 15 km s^ 1 . (b) NGC 3184. Contour levels are 5, 10, 15 km s^ 1 and the colour scale peak 
is 25 kms" 1 . (c) NGC 3938. Contour levels and colour scale are the same as for NGC 3184. (d) NGC 4736. Contour levels are 10, 20, 
... 60 km s _1 and the colour scale peak is 85 km s _1 . (e) NGC 4826. Contour levels and colour scale are the same as for NGC 4736. (f) 
NGC 5055. Contour levels and colour scale are the same as for NGC 4736. 
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Figure 2. CO J=3-2 velocity dispersion for five galaxies which are members of the Virgo cluster and NGC 2403, which is the closest 
galaxy in our sample. Colour scale runs from (black) to 50 km s _1 (white) and contour levels are 10, 20, 30, 40, 50 km s" 1 unless 
otherwise noted, (a) NGC 4254. (b) NGC 4303. (c) NGC 4321. (d) NGC 4501. Contour levels are 10, 20, ... 60 km s" 1 and the colour 
scale peak is 65 km s . (e) NGC 4535. (f) NGC 2403. Contour levels are 5, 10, 15 km s — 1 and the colour scale peak is 25 km s . 
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Table 1. Galaxy Properties and Velocity Dispersions 



Galaxy 


Vhel 


Type" 


D 25 a 


i a 


PA b 


K(tot) c 


log L FIR d 


Distance 6 




° c — c 


°Hi h 




(km s — ) 




(') 


(°) 


(°) 


(mag) 


(L©) 


(Mpc) 


(km s _1 ) 


(km s — 1 ) 


(km s ) 


NGC 628 


648 


SA(s)c 


10.5 


7 


20 


6.84 


9.55 


7.3 


4.1 ± 2.5 


3.5 ± 0.1 


12 


NGC 2403 


121 


SAB(s)cd 


21.9 


56 


124 


6.19 


9.10 


3.2 


5.2 ±4.9 


3.9 ± 0.2 


14 


NGC 3184 


574 


SAB(rs)cd 


7.4 


16 


179 


7.78 


9.69 


11.1 


6.8 ± 5.9 


5.8 ± 0.4 


13 


NGC 3938 


829 


SA(s)c 


5.4 


25 


52 


7.81 


9.72 


14.2 


4.4 ± 2.9 


2.7 ±0.3 


10' 


NGC 4254 


2412 


Sa(s)c 


5.4 


29 


56 


6.93 


10.62 


16.7 


9.2 ±4.5 


8.5 ± 0.3 




NGC 4303 


1577 


SAB(rs)bc 


6.5 


17 


318 


6.84 


10.59 


16.7 


8.0 ±4.2 


7.2 ± 0.2 




NGC 4321 


1599 


SAB(s)bc 


7.6 


32 


143 


6.59 


10.47 


16.7 


12 ± 11 


11.1 ±0.6 




NGC 4501 


2213 


SA(rs)b 


6.9 


58 


140 


6.27 


10.41 


16.7 


20 ± 15 


20 ±2 




NGC 4535 


1972 


SAB(s)c 


7.1 


41 





7.38 


10.12 


16.7 


5.1 ± 4.4 


3.7 ±0.4 




NGC 4736 


295 


(R)SA(r)ab 


11.2 


36 


296 


5.11 


9.79 


5.2 


12 ±9 


11. 9± 0.7 


25 


NGC 4826 


390 


(R)SA(rs)ab 


10.0 


57 


121 


5.33 


9.85 


7.5 


14 ± 11 


13 ±2 


60 


NGC 5055 


495 


SA(rs)bc 


12.6 


59 


102 


5.61 


10.30 


8.0 


9.4 ± 6.6 


8.7 ±0.3 


24 



a Buta et al. (2007) 

b Position angle used for radial averages. NGC 628, NGC 3184: Tamburro et al. (2008) : NGC 2403, NGC 4736, NGC 4826, 

NGC 5055: de Blok et al. (2008) ; NGC 3938: Paturel et al. (2000) : NGC 4254, NGC 4303, NGC 4501, NGC 4535: 

(Cavatte et al. 1990) : NGC 4321: Knapen et al. (1993) . 

c K-band total apparent magnitude from Jarrett et al. (2003). 

d From Sanders et al. (2003) adjusted for the distances adopted here. 

6 NGC 628: Karachentsev et al. (2004) : NGC 2403: Freedman et al. (2001) : NGC3184: Leonard et al. (2002) : NGC 3938, 
NGC 5055: distance calculated using Hubble flow distance with velocity corrected for Virgo infall (Mould et al. 2000) and 
H = 73 km s" 1 Mpc" 1 ; Virgo Cluster: Mei et al. (2007) : NGC 4736, NGC 4826: Tonrv et al. (2001) . 

* Observed CO J=3-2 velocity dispersion and standard deviation. Regions of high velocity dispersion in the outer portions 

of the disk were masked before calculating the velocity dispersion (see text). Excludes a 5.3 kpc diameter region in the centre of each 

Virgo galaxy (9 pixels), NGC 3938 (11 pixels) and NGC 5055 (19 pixels). All pixels are included for NGC 628, NGC2403, 

and NGC 3184. For NGC 4736 and NGC 4826, a 9 pixel region is excluded to give an upper limit to the outer disk velocity dispersion 

(see text). Note that the uncertainty in the mean for a ooa is the same as for <r c _ c . 

9 Mean cloud-cloud velocity dispersion and the uncertainty in the mean corrected for the contribution from the internal 
velocity dispersion of individual giant molecular clouds (see text). Note that the uncertainty in the mean is calculated from the 
standard deviation given in the previous column by dividing by vJV, where N is the number of measurements for a given galaxy. 
h HI velocity dispersion measured from naturally weighted maps from Waltej^et^iL^ (2008) over same area as the CO 
measurement. 

'From van der Kruit fc Shostak (1982) . 



sample. Beam smearing is expected to have an effect only 
on the measured velocity dispersions for the central pixel in 
each plot. Three of the galaxies (NGC 628, NGC 2403, and 
NGC 3184) show very flat profiles as a function of radius. 
These late type spirals are also the least luminous and so 
presumably are also the least massive galaxies in our sam- 
ple. The remaining 9 galaxies show a central peak in the 
velocity dispersion extending to 0.2 — 0.4r25. We attribute 
part of this increase to the effects of a steeply rising rotation 
curve at our relatively low (0.2-1.2 kpc) spatial resolution in 
these galaxies. Thus, for these galaxies, we measure the aver- 
age velocity dispersion in the outer disk excluding a central 
circular region. The exclusion aperture was chosen by ex- 
amining the individual velocity dispersion maps to identify 
the smallest aperture that would exclude the region of most 
obviously enhanced dispersion. We found that a diameter of 
5.3 kpc (9 pixels or 65" at the distance of Virgo) worked well 
for all five galaxies in the Virgo cluster and that a similar 
physical region was appropriate for the other galaxies in our 
sample. However, the maps of NGC 4736 and NGC 4826 
contain no data at larger radii, and so for these two galaxies 
we quote disk velocity dispersions excluding only a central 9 
pixel diameter region. The disk velocity dispersions of these 
two galaxies are likely biased to be somewhat higher than 
those of the other galaxies in our sample and thus we do not 



tween our sample and the THINGS sample l|Walter et al.l 
2008), we can compare the velocity dispersions in the atomic 
and molecular gas directly. The moment 2 maps for the 
THINGS sample we re produced using a slightly different 
masking technique ( Walter et alj 120081 ) than the one we 
adopted for our CO analysis. In the THINGS processing, the 
data are first smoothed to a resolution of 30" . A mask is then 
made keeping only those pixels and velocity channels (width 
either 2.6 or 5.2 km s _1 ) where the emission exceeds 2er in 
3 adjacent velocity channels. We used the natural weighted 
maps of the HI velocity dispersion and measure the disper- 
sion over the same region in which the CO dispersions have 
been measured. These velocity dispersions are also given in 
Table [T] Note that these HI dispersions are som ewhat larger 
than t he typical value of 11 ± 3 km s _1 given in lLerov et al.l 
(2008) because they are measured in the inner rather than 
the outer disks of the galaxies. 



3 VELOCITY DISPERSIONS IN THE 
MOLECULAR GAS DISK 

The observed CO velocity dispersion maps are shown in Fig- 
ures Q] and [2] Figures [3] and [4] show the observed velocity 
dispersion as a function of radius for the 12 galaxies in our 
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R A 25 

Figure 3. Observed CO J=3-2 velocity dispersion as a function of radius normalized by r25. These velocity dispersions have not been 
corrected for any contribution from the internal velocity dispersion of individual giant molecular clouds (see text). Where appropriate, 
the central region excluded from calculating the disk-avcraged velocity dispersion is indicated by the dashed line. 

I 1 I 1 I 1 I 1 
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R/r 25 

Figure 4. Observed CO J =3-2 velocity dispersion as a function of radius normalized by r25, showing nine galaxies on a single plot. 
Data points with r < 2.65 kpc are not plotted for the 6 galaxies with central peaks in the velocity dispersion. 
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include them in calculating the average outer disk velocity 
dispersion. 

The average observed velocity dispersions in the molec- 
ular gas measured for the 12 galaxies in our sample are given 
in Table Q] These velocity dispersions are consistent with 
previ ous meas urements for the Milky W ay l|Stark k, Brand! 
Il989h . M 33, JWilson fe Scoville! Il990l ) and other nearb y 
galaxies l|Combes fc Becauaerj Il997l ; IWalsh et al l |2002| ). 
The observed velocity dispersion in the molecular gas disk 
ranges from a low of 4.1 km s _1 in NGC 628 to a high of 
20.1 km s" 1 in NGC 4501. We exclude NGC 4736 and NGC 
4826 from our averages because we can only trace the veloc- 
ity dispersions in the inner region where the profile is still 
rising steeply. We further exclude NGC 4501 because its very 
high velocity dispersion and evidence for double peaked line 
profiles in the north-west portion of the disk suggest that 
the molecular gas has b een affected by the s ame ram pres- 
sure effects seen in HI (|Vollmer et al.|[2008l ). The average 
observed value for the remaining 9 galaxies is 7.1 ± 0.9 km 
s _1 (standard deviation 2.6 km s _1 ), 50% smaller than the 
average of 11 ± 3 km s _1 measur ed for the atomic gas in 
the outer disks of spiral galaxies (|Lerov et al.ll2008h . If we 
compare the atomic and molecular velocity dispersions in 
the same region of the disk (Table [TJ, the observed veloc- 
ity dispersion in the molecular gas is on average ~ 2 times 
smaller than the dispersion in the atomic gas. 



3.1 The effect of cloud internal velocity 
dispersions 

The low values measured for the velocity dispersion, par- 
ticularly for NGC 628 and NGC 3938, lead us to consider 
the effect of internal velocity dispersions of individual giant 
molecular clouds on the observed velocity dispersion. A sec- 
ond complicating effect, that of an anisotropy between the 
vertical and in-plane motions of the GMCs, would be ex- 
pected to produce a correlation of velocity dispersion with 
inclination, which is not seen in our data (Figure [Bl]). Nev- 
ertheless, we discuss the potential magnitude of this effect 
in Appendix B. 

Individual giant molecular clouds have internal velocity 
dispersions that are substantially larger than the thermal 
line-widths expected for gas with physic a l tem peratures of 
10-30 K. In the study of ISolomon et al.l (| 19871 1. individual 
GMCs show internal velocity dispersions ranging from 1- 
8 km s _ . These internal velocity dispersions are obtained 
from an intensity-weighted measure of the dispersion in ra- 
dial velocity within a single cloud and thus are directly com- 
parable (in technique) to the intensity-weighted dispersions 
averaged over the galactic disk presented here. What we 
need is an estimate of the contribution of the internal ve- 
locity dispersions of an ensemble of GMCs with a range of 
masses to the observed velocity dispersion in the disk. With 
this value, it is easy to show that the observed velocity dis- 
persion, a v , is given by o\, — Oc-c + a\ nt where aint is the 
mass-weighted average internal velocity dispersion for our 
ensemble of clouds and <t c _ c is the value we are interested 
in measuring, that is, the cloud-cloud velocity dispersion. 

It is well known that the cloud internal velocity dis- 
persion, ai n t, rel ates to the cloud mass as M oc af nt , as 
first proposed bvlHenriksen fc Turner] ll 19841 '). Observation 



dispersion aint (in km s ) relates to the cloud mass M (in 
M Q ) as M = 2000cri t . If we assume a cloud mass function 
dN/dm oc m~ a , then we can estimate the mass- weighted av- 
erage internal velocity dispersion for an ensemble of clouds 
between Mi om and Mhigh- Since colu mn density is prop or- 
tional to intensity for the CO lines (|Strong et al.|[l98g ). a 
mass-weighted average is also equivalent to an intensity- 
weighted average. The mass-weighted square of the internal 
velocity dispersion is then J ma^ nt dN/ J md N which, using 
the re lation between a^t and M given by ISolomon et al.l 
(|l987t ) and assuming a/2, reduces to 



lvl high 



low 



44.7 2.5 - a 



M 2 ~ a 

lvl high 



Mf. 



For a = 1.5 [|Sanders et ail Il985h. M hiah = 10 6 M and 
M iow = 10 4 M Q (e.g., ISolomon etal"1ll987T ). a in t = 3.5 km 



hrst proposed by ricnrikscn fc turner! l|19o41 1. Ubservation- 
ally, (Solomon et alj ( 1987 ) find that the internal velocity 



s~ . The exact result clearly depends on the details of the 
model. 

Interestingly, this estimate of the average cloud internal 
velocity dispersion is comparable to the value of 4.1 km s" 1 
observed in the disk of NGC 628. The average CO J=3-2 
intensity in NGC 628 of 0.8 K km s _1 is the lowest inten- 
sity in our sample and translates into an average mass per 
JCMT beam of only 4 x 10 6 M Q . For a GMC mass func- 
tion with slope a — 1.5, this surface density translates into 
an average of just 10 clouds with 10 5 < M < 10 6 M per 
beam. This calculation suggests that, although the number 
of massive (> 10 Mq) GMCs per beam may be subject to 
small number statistics in NGC 628, their effect on the veloc- 
ity dispersion should be observable on average if such large 
clouds are present in NGC 628. This analysis suggests that 
NGC 628 is likely deficient the largest molecular clouds; if 
we instead limit the maximum cloud mass to 10 Mq, then 
the mass-weighted average internal velocity dispersion re- 
duces to a^t = 2.2 km s _1 . The larger average integrated 
intensities and observed velocity dispersions suggests that 
the molecular cloud mass spectrum is likely to be fully pop- 
ulated in the remaining galaxies in our sample. Thus, in 
correcting the observed velocity dispersions for the effects 
of internal cloud velocity dispersions, we use a value for aint 
of 2.2 km s _1 for NGC 628 and 3.5 km s _1 for the rest of 
the sample. The corrected values for the cloud-cloud velocity 
dispersion are given in Table [1] along with the uncertainty 
in the mean cloud-cloud velocity dispersion for each galaxy. 
The uncertainty in the mean has been calculated assuming 
that the errors are normally distributed, so that the uncer- 
tainty in the mean is given by the standard deviation divided 
by \/~N , where N is the number of measurements for a given 
galaxy. 



3.2 The true vertical velocity dispersion of the 
molecular gas 

Based on the discussion in the previous section, the ob- 
served velocity dispersion includes both a contribution from 
the internal velocity dispersion of individual giant molecular 
clouds as well as a contribution from the cloud-cloud velocity 
dispersion. It is this second quantity, the cloud-cloud veloc- 
ity dispersion, which is of interest for determining the molec- 
ular gas scale height, analysing the stability of the gas disk, 
etc. The small observed velocity dispersions imply that the 
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correction for the internal velocity dispersion of the clouds 
is not negligible. Correcting for the average internal cloud 
velocity dispersion as discussed in £13.11 gives measures of 
the cloud-cloud velocity dispersions which range from 2.7 to 
19.8 km s _1 . Again excluding NGC 4501, NGC 4736, and 
NGC 4826, we obtain an average value of 6.1 ± 1.0 km s -1 
(standard deviation 2.9 km s _1 ). 



4 DISCUSSION 

4.1 Molecular and atomic gas velocity dispersions 

The average cloud-cloud velocity dispersion in the molecu- 
lar gas in this group of spiral galaxies of 6.1 ± 1.0 km s _1 is 
slightly smaller than previous estimates for external galaxies 
which di d not take the effect of internal cloud dispersion into 
account (|Combes fc Becauaertll 19971 ; I Walsh et al . 2002), but 
is comparable to recent measurements for th e cloud-cloud 
dispe rsion of mass ive clouds in the Gala xy l|Stark fc Led 
120061 ) and in M33 (|Wilson fc Scovillelll990fl . 

4-1.1 Effects of CO optical depth 

An additional possible concern in interpreting the CO data is 
the high optical depth of the CO line. Since we are interested 
in the cloud-cloud velocity dispersion, it is important to 
check whether the surface density of the molecular gas is suf- 
ficiently high that shielding of one cloud by another at a sim- 
ilar velocity could affect our measurements. For all but one 
of the galaxies in our sample, the projected surface density 
Nh is substantially smaller than the typical surface density 
of indi vidual giant m olecular clouds (Nh — 1.5 ± 0.3 x 10 22 
cm -\ lMcKeelll989h . and thus cloud-cloud shielding is un- 
likely to be a problem. Only in NGC 4736, where the surface 
density may be as high as 3 x 10 22 cm -2 (assuming a CO 
J=3-2/J=l-0 line ratio of 0.3 and a CO to H2 convers ion 
factor of 2 x 10 20 cm -2 (K km s" 1 )" 1 IStrong et af]| 19881 ), is 
it possible the clouds may shield one another to some degree. 
However, the actual surface density depends sensitively on 
the value that we assume for the CO J=3-2/J=l-0 line ra- 
tio; if we adopt a value of 0.6 which is more appropriate for 
the dense gas, rather than the value of 0.3 found in a direct 
comparison of CO J=3-2 and J=l-0 ma ps of NGLS galaxies 
l|Wilson et al.ll2009l ; [Warren et al.ll2010f ), then the projected 
disk-averaged column density is similar to that of individual 
GMCs. Thus, we conclude that we can ignore the possibility 
of cloud-cloud shielding in the disks of the galaxies in this 
sample. 

4-1.2 Observed HI velocity dispersions 

ICombes fc BecquaertJ |l997h showed that the molecular and 
atomic components in NGC 628 and NGC 3938 had verti- 
cal velocity dispersions that are similar to each other and 
virtually constant with radius. As a consequence, the au- 
thors suggested that the two layers had similar scale heights 
and that HI and H2 simply represent different phases of the 
same dynamical component. It is clear, however, from ob- 
servations of our own Milky Way as well as edge-on galaxies 
that the molecular disk is significantly thinner than the HI 
disk. Our observed velocity dispersions are comparable to 



those found by ICombes fc Becauaerj l|l997f ) , although we 
have shown that it is necessary to correct for internal ve- 
locity dispersion to determine the cloud-cloud velocity dis- 
persion that is relevant for scale height determinations. For 
HI, on the other hand, we know that its stable phases occur 
in both hot diffuse gas as well as colder clouds (the warm 
neutral medium and cold neutral medium, respectively). Al- 
though discrete HI clouds are id entifiable in our Galaxy (see 
iLockwoodlkoOltlStil et al.ll2006T ). the clouds are smaller and 
warmer than molecular clouds. Because of the admixture of 
cool and warm HI components over the size scale of the beam 
in external galaxies, similar corrections need not be applied 
to the HI velocity dispersions. Also, since the HI is a diffuse 
collisional me dium, its veloci ty dispersion is assumed to be 
isotropic fe.g. lMalhotra|[T995l ). We have therefore compared 
the observed HI dispersion directly with our cloud-cloud ve- 
locity dispersion derived from the CO data. 

4-1.3 Comparing CO and HI velocity dispersions 

The published HI v elocity dispersion map for NGC 4501 
l|Vollmer et al.ll2008l ) shows velocity dispersions in the range 
of 10-20 km s _1 across the extended disk, comparable to the 
cloud-cloud velocity dispersion (Table[T} over the same area. 
High resolution HI d ata have been p ublished for NGC 432 1 
l|Knapen et al.lll993l ) and NGC 4254 l|Phookun et al.lll993T l. 
but the velocity dispersion maps are not presented; we could 
find no published high resolution HI data for NGC 4303 and 
NGC 4535. An older measurement of the HI velocity disper- 
sion of 10 km s _1 for NGC 3938 l|van der Kruit fc Shostakl 
1982) is roughly 3 times larger than our derived cloud-cloud 
velocity dispersion. 

We have compared the velocity dispersions in the molec- 
ular and atomic gas measured over the same regions for six 
galaxies in our sample for which maps of the atomic gas 
velocity dispersion are available from the THINGS sample 
l| Walter et al.l l2008). On average, the cloud-cloud velocity 
dispersion determined from CO data is roughly half that of 
the atomic gas. The exception to this trend is NGC 4826, 
for which the cloud-cloud velocity dispersion is a factor of 
5 times smaller than the the HI velocity dispersion. How- 
ever, this galaxy has quite a steep radial gradient in the ob- 
served velocity dispersion (Figure [3} and so the average disk 
value depends sensitively on what aperture is used. Overall, 
these results imply that the scale height of the molecular 
gas is roughly half that of the atomic gas. This estimate 
of the ratio of the scale heights is roughly consistent with 
the estimate of the relativ e HI and CO scale heights i n the 
Milky Way as discussed in lCombes fc BecquaertJ (|l997l ). We 
should note, however, that narrow (4-5 km s _1 ) HI line com- 
ponents have been seen in a few Local Group dwarf galax- 
ies observed with high spatial resolution (lYoung fc Lolll996l . 
1 19971 ; Ide Blok fc Walter! liooe} ) . These narrow lines, which 
may represent atomic gas cooling to form the precursors of 
GMCs, are typically superimposed on a broader underlying 
line which likely represents the overall warm neutral atomic 
medium. 

4.2 Correlations with other global properties 

The infrared luminosity is a us eful tracer of the star forma- 
tion activity within a galaxy l|Kennicut j 1 1998h ■ Enhanced 
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star formation might be expected to increase the veloc- 
ity dispersion in the interstellar m edium via stellar winds 
and supernova explosions. Indeed, iTamburro et~ai1 (2009) 
find evidence that supernovae linked to recent star forma- 
tion are important for maintaining the velocity dispersion 
in the atomic gas. We find a statistically significant cor- 
relation (95% confidence level) of both the observed and 
the cloud-cloud velo city dispersion with infrared luminosity 
l|Sanders et al.| [2003') for the 9 galaxies in our sample. Thus, 
like the atomic gas, there is some evidence for star forma- 
tion activity enhancing the velocity dispersion in the dense 
molecular gas. We also find a statistically significant corre- 
lation (95% confidence level) of both the observed and the 
cloud-cloud velocity dispersion w ith the absolute K -band 
and absolute B-band mag nitude (jjarrett et alJl2003h . As- 
suming the K-band magnitude traces the total stellar mass, 
this correlation suggests that the cloud-cloud velocity dis- 
persion is enhanced in more masssive galaxies. For all these 
correlations, however, it is important to keep in mind that 
we have compared the average velocity dispersion in the 
disk (excluding the central regions) with a global luminos- 
ity which would include the contribution from the central 
regions. 

There also appears to be a dependence in our sample 
of the velocity dispersion on galaxy morphology. We divided 
our sample into early (type ab-bc) and late (type c-cd) spi- 
rals and calculated the average velocity dispersion for each 
sub-sample. The average observed dispersion for 6 late-type 
spirals is 5.8 ± 0.3 km s~\ while the value for 3 early- type 
spirals is 9.7 ± 0.6 km s" 1 . These two values differ at the 6 
sigma level. Of course, this dependence on morphology could 
be tracing some other variable in our sample, such as mass 
or star formation activity. In fact, the average infrared lu- 
minosities differ by a factor of three between the early- and 
late-type samples, while the K-band luminosity differs by a 
factor of four. Thus, it seems most likely that the depen- 
dence on morphology is driven by differences in the average 
star formation rates and masses between the early and late 
type samples. 



4.3 Disk stability and star formation laws 

The velocity dispersion of the interstellar medium is an 
important parameter in understanding t he stability of 
alactic disks against gravitational collapse. I Jog fc Solomonl 
1984al lbl) showed that a rotating disk composed of two flu- 
ids with different surface densities and velocity dispersions 
can be unstable to perturbations even if each fluid disk is 
stable on its own. They also showed that a galaxy must be 
treated as a two fluid system when the colder gaseous com- 
ponent comprises as little as 10% of the total mass of the 
system. Treating the stellar component as a collisional fluid, 
they found that the critical velocity dispersion required for 
stability in this higher dispersion system is larger in the 
presence of a secon d cold fluid tha n if the stellar system is 
treated in isolation. Rafikov (2001) considered the case of a 
single fluid and multiple collisionless (stellar) components, 
but found that the stability conditions were quite similar in 
the two fluid case whether or not the stellar component was 

trea ted as collis io nless. 

| Yang et al.l (|2007l ) applied the approach of Rafikov 
(|200lF to compare the location of star forming regions in the 



Large Magellanic Cloud (LMC) with the disk stability cri- 
terion. They adopted a gas velocity dispersion of 5 km s~\ 
very similar to the value found in our sample of spiral galax- 
ies. They found that there was a much better spatial correla- 
tion of star forming regions with regions where the disk was 
unstable when instability was ev aluated using the c ontribu- 
tion of both the gas and the stars. iLerov et al. (l2008h applied 



the two collisional fluid equation of lRafikovl ll200ll) to a num - 
ber of galaxies in the THINGS sample (| Walter et al.|[2008h . 
They adopted a gas velocity dispersion of 11 k m s -1 from 
measu rements of the HI component. Similarly to lYang et al.l 
(2007), they found that the Q values were smaller when both 
the stars and the gas were included in the analysis; however, 
they found that the disks were globally stable against large 
scale perturbations. 



The results for the molecular velocity dispersion pre- 
sented here imply that the interstellar medium in spiral 
galaxies is itself a multi-fluid system, with the molecular gas 
disk being dynamically significantly colder than the atomic 
gas disk. Since previous analyses have found that the colder 
of t he two components has the larger effect on the disk stabil- 
ity (|Jog fc Solomonlll984ar ). it seems likely that, if we were to 
use one phase of the ISM in our stability analysis, it should 
be the dynamically coldest phase, which is the molecular 
gas. This would be especially true for the inner regions of 
those spiral disks where the ISM mass is predominantly H2. 
The Toomre Q parameter for the combined system tends to 
decrease as the velocity dispersion of the coldest (gas) com- 
ponent d ecreases and as the ga s mass fractio n (E ga s/£*) 
increases (|jog fc Solomonlll984bl ; lRafikovll200ll). This e ffect 
suggests that the Q values derived bv lLerov et al.l |2008T l are 
likely to be significantly larger than the values that would be 
derived if the molecular gas values for surface density and 
velocity dispersion we re used in t h e ana lysis, a possibility 
that was discussed in ILerov et al.l (|2008r i as well. It would 
be informative to repeat the stability analysis for the galax- 
ies in common between our sample and the THINGS sample 
to see if the disks remain globally stable as judged by the Q 
parameter when the cloud-cloud velocity dispersion of the 
molecular gas is included. We plan to investigate this issue 
further in a future paper. 



The need to treat the ISM itself as a multi-fluid system 
suggests that we may not be able to separate the discussion 
of the instabilities in the atomic and molecular phases. For 
example, we might imagine that the Q value for the atomic 
gas governs the formation of molecular clouds from the HI 
disk, while the star formation from those clouds would then 
be governed by the Q value f or the molecula r disk. How- 
ever, the two fluid analysis by I Jog fc Solomonl (|l984aT ) and 
I Rafikov! (|2001i ) implies that the Q value and hence the stabil- 
ity properties of the atomic disk are affected by the presence 
of the dynamically colder molecular disk. In particular, an 
atomic gas disk with a stable value of Q when considered in 
isolation might be found to be unstable when the presence 
of a relatively small quantity of molecular gas is included. Of 
course, a complete analysis would include the prop erties of 
the s tellar disk (perhaps multiple components, as in [Rafikov 
2001) in addition to the two fluid components of the ISM. 
Such an analysis is beyond the scope of this paper. 



The JCMT Nearby Galaxies Legacy Survey IV. Velocity Dispersions in the Molecular Gas 11 



4.4 Implications for turbulent support of the ISM 

The scale height of a gaseous component in hydrostatic equi- 
librium depends on both gravitational and pressure gradi- 
ents. Gravitational gradients originate in the disk as well as 
a dark matter halo potential. The pressure will include inter- 
nal thermal pressure, turbulent pressure and, if important, 
pressure due to cosmic rays and magnetic fields. At the ra- 
dius of the Sun in the Milky Way, the magnetic and cosmic 
ray pressures supply about half of the total pressure and the 
thermal and turbulent pressures supply the remaining half 
(Hanasz & Lesch 2004). At the temperatures of the molecu- 
lar components, turbulent pressure domin ates over the rmal 
pressure within and between clouds (e.g. iBruntl 120031 ). In 
the HI, thermal pressure likely dominates in the Warm Neu- 
tral Medium whe reas both may be impo rtant for the Cold 
Neutral Medium (jHeiles fc Trolandl [20031 '). Measured veloc- 
ity dispersions can only probe the dynamical components 
(thermal and turbulent) and therefore can only predict the 
scale height of the gas provided that the dynamical compo- 
nents dominate over the magnetic and cosmic ray pressures 
and provided that the disk ISM is in pressure equilibrium. 

Star formation activity is typically associated with heat- 
ing sources and sources of turbulence through supernovae 
and stellar winds. Most of the galaxies in our sample show 
a radial decline in the CO velocity dispersion (Figure Q. 
The veloci ty dispersion in HI is also a declining function 
of radius (|Tamburro et all l2009h . although HI disks re- 
tain signi ficant velocity dispersions outside any star formin g 
disk (e.g. ISellwood fc Balbusl Il999l ; iTamburro et~ai]|2009t ). 
iTamburro et all (|2009h found a good correlation between the 
kinetic energy of the HI and the star formation rate and sug- 
gest that supernovae are likely sufficient to maintain the HI 
velocity dispersion in regions with significant star formation. 
The correlation of the cloud-cloud velocity dispersion with 
the star formation rate as traced by the infrared luminos- 
ity suggests that star formation may also be the dominant 
source of turbulence in the H2-rich parts of spiral disks. 

However, some galaxies do show a lack of correla- 
tion between turbulence in the atomi c gas and star for- 
mation as a function of radius, (e.g. Dickey et alj Il990l ; 
Ivan Zee fc Brvantl 1 19991 ; iPetric fc Rupenl l2007h ■ The three 
lowest luminosity galaxies in our sample show a roughly 
constant cloud- cloud velocity dispersion as a function of 
radius (see also ICombes fc Becquaertl Il997l ) These results 
suggest that non-thermal energy input in the form of turbu- 
lence that is unrelated to star formation may be important 
for understanding the vertical velocity dispersions in disks. 
Magnetic instabilities may be the most promis ing source 
of turbulence. Parker instabilities (|Parkerl Il966h are well- 
known, but seem to require c osmic rays from supern ovae as 
triggers to be most effective jHanasz fc Leschll2000l) . Thus, 
Parker instabilities cannot account for the observed disper- 
sions seen outside of the active star forming disk, but could 
be important for the inner molecular disks where star forma- 
tion is occurring. Another possibili ty is the Balbus-Hawley 
(or m agneto-rotational) instability (|Balbus fc Hawlevlll99ll . 



1998) , or other instabilitie s related to magne t ic stre sses (e.g 



Sellwood fc Balbusl [l999l ). ITamburro etail (2009) suggest 



instability by IPiontek fc Ostrikerl (|2007h for the two-phase 
HI component have shown that the turbulent velocity ampli- 
tude varies as oc n~ ' 77 , where n is the mean density of the 
component. If this turbulence is dominant in both molecu- 
lar and atomic gas, then the difference in velocity dispersion 
between the two components may simply relate to their re- 
spective densities. 

4.5 Disk galaxies at high redshift 

Observations of the molecular gas content in galaxies at 
redshifts 2 = 1 — 3 are often made using the CO J=3-2 
line, which at these redshifts falls into the standard millime- 
tre observing windows. Most high-resolution detections of 
high redshift galaxies are u nusually luminous sy stems, such 
as submillimetre ga laxies (|Tacconi et alj 120081 ) or quasars 
|Coppin et al.l l2008). Although many of these high redshift 
detections appear as simple point sources at current reso- 
lution limits, there are a few systems for which we can be- 
gin to resolve the velocity structure. Even more interesting 
from the point of view of this paper are the observations of 
more normal galaxies at redshifts 1-1.5, which are becom- 



ing ava ilable ( Daddi et alj 12008 



iDaddi et al]|2010l ; iTacconi et allFoiOh . These gal. 



Dannerbauer et all 2009; 



lax ics arc 



that a combination of thermal broadening and magneto- 
rotational instabilities can account for the HI velocity dis- 
persion beyond r25. Simulations of the magneto- rotational 



found to be quite extended (~ 10 kpc), signif icantly more 
so th a n the quasars and su bmillimeter galaxies l|Daddi et al.1 
|2010| ; ITacconi et al. I l201Ch Given the sensitivity and wide 
field of view of our observations, it is interesting to examine 
how our galaxies would appear if placed at a redshift of 1. 

We have processed a subset of our galaxies to see what 
changes are produced in an analysis of the velocity disper- 
sion when the spatial resolution of the data is degraded. For 
three of the brightest galaxies in the Virgo Cluster (NGC 
4254, NGC 4321, and NGC 4303), we convolved the baseline- 
subtracted data cube with a 54" gaussian to achieve an ef- 
fective beam of 56" or 4.5 kpc at a distance of 16.7 Mpc. 
We then constructed moment maps from the smoothed data 
cube using the methods described in §[2] A comparison of the 
resulting velocity dispersion maps for NGC 4254 is shown 
in Figure [S] Depending on how much of the central high 
dispersion regions is excluded from the analysis, the average 
velocity dispersion in the outer disk of the low resolution 
image is on average twice as large as the value in the high 
resolution image. Similar increases of about a factor of two 
in the average velocity dispersion in the disk are seen for 
NGC 4321 and NGC 4303. The images also suggest that a 
plot of velocity dispersion as a function of radius would yield 
higher values at a given radius in the low resolution image. 
These effects of resolution will need to be taken into account 
as we begin to accumulate data on the dense molecular gas 
properties in galaxies at high redshift. 

The average values of 20-30 km s~ seen at low res- 
olution in the disks of these three galaxies are quite simi- 
lar to the value of ~ 20 km s' 1 found in EGS13035123 by 
ITacconi et ail (|201Ch . EGS13035123 has 10-20 times the star 
formation rate, H2 gas mass, and stellar mass compared 
to NGC 4254 ilKennicutt et all 120031 ; IWilson et all 120091 : 
iKranz et aLll2003h . However. EGS13035123 is also roughly 
a factor of three la rger in its CO radius than NGC 4254 
l|Wilson et al.ll2009h , which implies that the gas surface den- 
sity, stellar surface density, and star formation rate surface 
density inside the molecular disks of these two galaxies are 
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Figure 5. (a) CO J=3-2 velocity dispersion for NGC 4254 at 
the native resolution of the JCMT data (14.5" or 1.15 kpc at a 
distance of 16.7 Mpc). Colour scale runs from to 60 km s _1 
and the contours are 10, 20, 30, 40, 50 km s _1 . (b) CO J=3-2 
velocity dispersion for NGC 4254 smoothed to a resolution of 4.5 
kpc (56") t o match the resolutio n of observations of disk galaxies 
at z ~ 1 of lTacconi et al" Colour scale and contours are 

the same as in (a). 



quite similar. Thus, the similar velocity dispersions between 
the two galaxies are likely related to the similar mass surface 
densities in the disk. The two galaxies also have quite sim- 
ilar molecular gas fractions of ~ 0.2 5 (|Tacconi et al.ll2010l ; 

1 Wilson et al.ll2009l ; lKranz et alj|2003h . This comparison sup- 
ports a picture where the high star formation rates seen at 

2 = 1 may be at least partly due to the presence of physically 
larger molecular gas disks at this epoch. 



5 CONCLUSIONS 

We have used large-area high velocity resolution CO J=3- 
2 observations of 12 nearby galaxies to study the verti- 
cal velocity dispersion in the dense molecular gas. Three 
of the galaxies show a roughly constant velocity disper- 
sion as a function of radius, while the other nine galaxies 
have a central peak followed by a fall-off with radius to 
typically 0.2 — 0.4>2s. Flat velocity dispersion profiles are 
seen only in some of the late-type spiral galaxies in our 
sample, and those with flat profiles have the lowest mass. 
The observed values of the velocity dispersion range from 
4.1 km s _1 to 20.1 km s _1 . These velocity dispersions are 
comparable to the internal velocity dispersions of individ- 
ual giant molecular clouds in our own and ot her galax- 
ies jSolomon et al.lll987l ; Iwilson fc Scovillelll99Ch . Correct- 
ing for these internal velocity dispersions yields an average 
cloud-cloud velocity dispersion of 6.1 ± 1.0 km s _1 mea- 
sured over 9 galaxies with good radial profiles. This cloud- 
cloud velocity dispersion is comparable to r ecent measure- 
ments in our own Galaxy (|Stark fc Lej2005l . l2006h and M33 
I Wilson fc ScoviUdll99ch . 

A direct comparison with the high resol ution HI maps 
from the THINGS survey | |Walter et alj|2008h in six galaxies 
show that the cloud-cloud velocity dispersion is on average 
twice as small as the velocity dispersion of the atomic gas, 
which implies a much smaller scale height for the molecular 
gas. Theoretical analyses of the stability of multi-component 
disks suggest that the dynamically coldes t component is 
the m os t important in driving instability (|jog fc SolomonI 
Il984al lbt iRafikovllioOll l . This analysis suggests that it is the 
properties of the dense molecular gas, rather than the atomic 
gas, that are the most important for determining whether 
galactic disks are stable against gravitational collapse, espe- 
cially where the mass of the ISM is H2 dominated. 

The cloud-cloud velocity dispersion is correlated at the 
95% confidence level with both the far-infrared luminosity 
and the K-band absolute magnitude. Thus, as for the atomic 
gas, we find evidence that star formation activity (as traced 
by the infrared luminosity) tends to increase the velocity 
dispersion in the dense molecular gas. The correlation with 
K magnitude, which traces the total stellar mass, implies 
that the cloud-cloud velocity dispersion is also enhanced in 
more massive galaxies. 

We have used our data to examine the apparent kine- 
matical properties of the molecular disk at a spatial reso- 
lution of 4.5 kpc chosen to match the best resolution for 
galaxies at redshifts 1-2 l|Tacconi et al.l[200ct l. A degrada- 
tion of the resolution from 1.2 kpc to 4.5 kpc results in an 
increase in the average velocity dispersion in the outer disk 
by a factor of two. The average velocity dispersion of NGC 
4254 viewed with 4.5 kpc resolution is quite comparable to 
the veloc ity dispersion of 20 km s _1 seen in a normal galaxy 
at z = 1 (jTamburro et al.ll2009l ). Both galaxies have compa- 
rable gas and stellar surface densities, as well as star forma- 
tion rate surface densities, which suggests that the higher 
star formation rates seen at z — 1 may be partly attributed 
to the presence of physically larger molecular disks. This 
analysis suggests that this data set can provide a valuable 
local bench mark in understanding lower spatial resolution 
observations of galaxies in the early universe. 
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APPENDIX A: COMPARISON WITH 
PREVIOUS PROCESSING METHODS 

There are two different methods that have been used to 
measure velocity dispersion s in molecular gas in galaxies. 
ICombes fc Becauaertl (|l997f ) measured velocity dispersions 
in NGC 628 and NGC 3 938 by fittin g gaus sian profiles to 
individual spectra, while IWalsh et al.l (2002) used moment 
2 maps of NGC 6946 to determine the average velocity dis- 
persion. While the moment 2 maps used in this paper are 
easier to compute and analyse than fitting individual profiles 
to many hundreds of spectra, they can potentially be sub- 
ject to some systematic effects depending on choices made in 
the processing. We investigate some of these effects here by 
comparing our data and moment 2 maps for NGC 628 and 
NGC 3938 with the results given by ICombes fc Becauaertl 
i| 19971 ) for the same galaxies. 

We first investigated whether the measured velocity dis- 
persion was affected by the angular and spectral resolution 



Table Al. Velocity dispersions measured with different spectral 
and spatial resolutions 



Data cube used" 


NGC 3938 


NGC 628 


Original 


3.7 


3.1 


Binned to 2.6 km s _1 


4.8 


3.9 


Convolved to 23" gaussian 


6.5 


4.2 


Convolved and binned 


7.4 


1.7 


Combes & Bccquacrt (1997) 


9 


6.5 



a All values measured for the CO J=3-2 line in km s . 

Velocity dispersion measured using a signal-to-noise cutoff of 3.5c. 



Table A2. Velocity dispersions measured with different signal- 
to-noise cutoffs for mask 



S/N cutoff used 


NGC 3938 


NGC 628 


5cr 


3.8 


3.4 


4.5o- 


4.8 


3.7 


4.0o- 


5.7 


3.9 


3.5o- 


6.5 


4.2 


3.O0- 


7.2 


4.6 


2.5ct 


9.1 


5.3 


Combes & Bccauacrt (1997) 


9 


6.5 



a All values measured for the CO J=3-2 line in km s . 



of the data. The CO J=l-0 data from lCombes fc Becauaertl 
l|l997l ) have an angular resolution of 23" and a spectral res- 
olution of 2.6 km s~\ while the JCMT CO J=3-2 data have 
an angular resolution of 14.5" and a spectral resolution of 
0.43 km s - . We produced three additional data cubes for 
each galaxy: one binned by 6 channels in velocity, one con- 
volved to achieve a gaussian beam of 23", and one that was 
both convolved and binned. From these cubes we made mo- 
ment 2 maps using the method described in §[2] but using a 
signal-to- noise cutoff for the mask of 3.5 sigma so as to in- 
clude only regions with relatively high signal to noise. The 
average velocity dispersion measured in each map is given in 
Table IA1I These results clearly show that data with poorer 
spectral or spatial resolution will give higher values for the 
velocity dispersion than data with better spectral or spa- 
tial resolution. The data also suggest that spatial resolu- 
tion may be a more important factor in increasing the mea- 
sured velocity dispersion than spectral resolution, as long 
as the spectral resolution is sufficient to resolve the lines. 
Interestingly, our average velocity dispersions for these two 
galaxies are only slig h tly sm aller than those measured by 
ICombes fc Beccmaert l|l997l ) when our data are smoothed 
spatially and spectrally to match the IRAM data. 

We next investigated the effect of the choice of the 
signal-to-noise cutoff for the masks used in creating the 
moment map. For this analysis, we used the data cubes 
that had been smoothed sp atially to match the data from 
ICombes fc Becauaertl (|l997l ) to give better signal-to-noise 
in our maps, but which had no spectral smooting applied. 
The results are given in Table lA"2l It is clear that the mea- 
sured velocity dispersion decreases systematically has higher 
signal-to-noise cutoffs are used in the mask. 

Finally, we performed gaussian fits to selected spectra 
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Table A3. Comparison of velocity dispersion from gaussian fits 
and second moment maps for selected positions in NGC 3938 



Position" 


Gaussian fit 


Moment 2 map 6 


(",") 


(km s _1 ) 


(km s —1 ) 


(0,44) 


9.4± 2.9 


7.8 


(0,22) 


14.0± 2.3 


19.3 


(0,0) 


23. Ozt 4.4 


19.3 


(0,-22) 


8.7± 1.3 


8.4 


(0,-44) 


15. 0± 3.2 


8.2 

















|;l | ft 




r i,i iftf — h 

1 yvf 1 






Iff 



a All values measured for the CO J=3-2 line. The spectra have been 
to a 23" beam and binned to 2.6 km s — 1 resolution. (0,0) correspond 
(11:52:49.6,44:07:17.7 J2000). 

b Measured from maps made using a 2.5ct cutoff. 

from the data cube of NGC 3938 that had been convolved 
to 23" and binned to 2.6 km s -1 resolution (Figure \A II and 
Table IA3I T o mimic approximatel y the selection of spec- 
tra shown in ICombes fe Becquaertl (|l997| ). we selected the 
spectrum with the highest velocity dispersion in a map made 
with a 3.5a cutoff, and then stepped away from this spec- 
trum in steps of 3 pixels (21.8") in the north and south 
directions. We obtain similar values for the velocity disper- 
sions from fits to the unbinned spectra; binned spectra are 
shown in Figure [AH for clarity. Although there is consider- 
able scatter in the velocity dispersions for individual spec- 
tra, the average value derived from gaussian fits is 1.2 times 
larger (the standard deviation 0.4) than the value derived 
from the moment 2 map. 

On the basis of this analysis, we opted to use a signal- 
to-noise cutoff of 2.5a on our original data cubes, as this 
choic e seemed to give good ag reement with the results 
from lCombes fc Becquaertl l|l997l ) when both data sets were 
matched in angular and frequency resolution. 



APPENDIX B: THE POSSIBLE EFFECT OF 
ANISOTROPIES IN THE VELOCITY 
DISPERSION 

There is no statistically significant correlation of velocity 
dispersion in the molecular gas with inclination in our data 
(Figure [Bill . The average observed velocity dispersion is 5.3 
km s _1 for the 4 galaxies with inclinations < 25° and 8.1 
km s _1 for the 8 galaxies with inclinations > 25° and these 
values agree within 2 sigma. We have also checked for any 
correlation of the velocity dispersion normalized by each of 
the star formation rate and D25 with inclination and again 
find no correlation. 

For an inclined galaxy, the observed velocity disper- 
sion is a combination of the in-plane (a r ,ae) and vertical 
velocity (a z ) dispersions. We will assume that ay = ag 
and will refer to the in-plane component of the velocity 
dispersion as a r from here on for simplicity. Thus, for a 
galaxy with inclination i, the observed velocity dispersion 



is a t s = V °* cos2 * + a r sm *• If t ne ra dial and vertical 
velocity dispersions of the molecular gas were isotropic, we 
would not expect any trend of velocity dispersion with incli- 
nation. If, however, the velocity ellipsoid of the giant molec- 
ular clouds (which contain most of the mass of the molec- 








1 


ik 






iff 1 Ti'[«||rn i|p j j 




Figure Al. Observed CO J=3-2 emission at five positions in 
NGC 3938. The spectra have been convolved to a 23" beam and 
binned to 2.6 km s — 1 resolution. Spectra are spaced by 22" along 
the declination axis ordered from north (top) to south (bottom) 
and the central spectrum corresponds to (11:52:49.6,44:07:17.7). 
Gaussian fits to each spectrum are overlaid; the bar indicates the 
region used to obtain the fit. 



ular ISM) is anisotropic, with a r > a z , then the observed 
velocity dispersion would tend to increase with increasing 
inclination. 

Unfortunately, there is no direct information on the 
shape of the velocity ellipsoid for giant molecular clouds. In 
theoretical models, the cloud veloc ity dispersion has been at- 
tribu ted to cloud-cloud scattering ( Ga mmie. Ostriker fc Jo j 
Il99ll) as well as to the driving effects of spiral structure 
( Thomasson et al.|[l99ll) . Since both these models were two- 
dimensional, they can give us no guidance on the relative 
strength of the vertical and in-plane velocity dispersions. 
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Figure Bl. Observed CO J=3-2 velocity dispersion as a func- 
tion of galaxy inclination angle. Error bars show the standard 
deviation of the observed values within an individual galaxy. 



In contrast, there ha ve been a number of studies of the 
stellar velocity ellipso id (Dclhave 1 19651 : iDehnen fc BinnevI 
ll998l ; lElias et alj|2006h ■ If the velocity ellipsoid of any stars 
can give us a clue to that of GMCs, it will be the youngest 
O and B type stars, which may be sufficiently young to still 
trace the motions of their parent clouds . The smallest v alues 
for the velocity ellipsoid can be found in |Pelhavell|l965h . who 
measured a relativ e velocity dispersion a z /< j r =0.75 for 35 
0-B5 supergiants. IDehnen fc BinnevI (1 19981 ) measure a sig- 
nificantly smaller value of <r z /a r = 0.38ig'io f° r 500 main 
sequence stars with —0.24 < B — V < 0.14. However, this 
colour range extends well into the A star range where main 
sequence lifetimes are larger than 100 Myr. Since this veloc- 
ity dispersion ratio is wel l known to decrease with the age 
of th e stellar population l|Delhavelll965l ; IDehnen fc BinnevI 
1998), this value may not be an appropriate one to adopt for 
GMCs. In a recen t study of OB stars using Hipparcos data, 
lElias et all <|2006l ) obtained a value of a z /a r = 0.64±0.06 for 
800 stars with spectral types 0-B6 within 1 kpc of the Sun. 
Dividing the sample into stars in the Gould's Belt and the 
Local Galactic Disk results in values of a z /a r = 0.72 ± 0.08 
and a z /o~r = 0-57 ± 0.09, respective ly. 

iGammie. Ostriker fc Jog |l99ll ) caution that a critical 
difference between star-cloud and cloud-cloud scattering, 
namely the relative sizes of the epicyclic amplitude and the 
cloud tidal radius, prevents direct application of stellar re- 
sults to molecular clouds. However, this model ignores the 
possible effect of spiral arms, which could conceivably intro- 
duce an anisotropy into the cloud motions. In estimating the 
possible correction to our observed velocity dispersions for 
the effects of anisotropy in the velocity ellipsoid, we adopt 

a conservative va lue of <J z lo r 0.6 , which is the 1<t upper 

limit obtained bv lElias et alj (2006) for their entire sample. 
With this value, corrections for the velocity ellipsoid only 
exceed 10% for inclinations greater than 27°. 



